home *** CD-ROM | disk | FTP | other *** search
/ Linux Cubed Series 3: Developer Tools / Linux Cubed Series 3 - Developer Tools.iso / devel / lang / fortran / libf2c-0.000 / libf2c-0 / libf2c / libF77 / pow_ri.c < prev    next >
Encoding:
C/C++ Source or Header  |  1994-03-20  |  348 b   |  35 lines
  1. #include "f2c.h"
  2.  
  3. #ifdef KR_headers
  4. double pow_ri(ap, bp) real *ap; integer *bp;
  5. #else
  6. double pow_ri(real *ap, integer *bp)
  7. #endif
  8. {
  9. double pow, x;
  10. integer n;
  11.  
  12. pow = 1;
  13. x = *ap;
  14. n = *bp;
  15.  
  16. if(n != 0)
  17.     {
  18.     if(n < 0)
  19.         {
  20.         n = -n;
  21.         x = 1/x;
  22.         }
  23.     for( ; ; )
  24.         {
  25.         if(n & 01)
  26.             pow *= x;
  27.         if(n >>= 1)
  28.             x *= x;
  29.         else
  30.             break;
  31.         }
  32.     }
  33. return(pow);
  34. }
  35.